Preprint typeset in JHEP style - HYPER VERSION 



Preheating After Modular Inflation 



Neil Barnaby 1 , J. Richard Bond 2 , Zhiqi Huang 3 and Lev Kofman 4 

Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 St. George 
St. Toronto, Ontario M5S 3H8 Canada 

1 barnaby@cita.utoronto .ca, 2 bond@cita.utor onto. ca, 3 zqhuang@astro.utoronto.ca, 
4 kofman@cita. utoronto. ca 



Abstract: We study (p)reheating in modular (closed string) inflationary scenarios, with a 
special emphasis on Kahler moduli/Roulette models. It is usually assumed that reheating 
in such models occurs through perturbative decays. However, we find that there are very 
strong non-perturbative preheating decay channels related to the particular shape of the 
inflaton potential (which is highly nonlinear and has a very steep minimum). Preheating 
after modular inflation, proceeding through a combination of tachyonic instability and 
broad-band parametric resonance, is perhaps the most violent example of preheating after 
inflation known in the literature. Further, we consider the subsequent transfer of energy to 
the standard model sector in scenarios where the standard model particles are confined to a 
D7-brane wrapping the inflationary blow-up cycle of the compactification manifold or, more 
interestingly, a non-inflationary blow-up cycle. We explicitly identify the decay channels 
of the inflaton in these two scenarios. We also consider the case where the inflationary 
cycle shrinks to the string scale at the end of inflation; here a field theoretical treatment of 
reheating is insufficient and one must turn instead to a stringy description. We estimate 
the decay rate of the inflaton and the reheat temperature for various scenarios. 
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1. Introduction 

Reheating at the endpoint of inflation is a crucial requirement for any successful model. 
Depending on how inflation ends and how the inflaton interacts, the process of reheating 
- the creation of particles from the decaying inflaton and subsequent thermalization - 
can proceed differently. Examples are known of perturbative inflaton reheating non- 
perturbative preheating [Q, |3| (leading to excitations of both bosons and fermions Q), 
and also reheating via string theory mechanisms such as intermediate Kaluza-Klein (KK) 
modes []|, |6|, etc. 

In the last several years there has been significant progress in string theory related 
to the realization of realistic compactifications with stabilized moduli |?], ^, ^, [Tcfl . This 
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progress has stimulated the development of a new generation of inflationary models based 
on such stabilized string theory constructions; see [11, |l^, [ll|] for reviews. Among the 
various possible string theory models of inflation, one can distinguish different classes 
depending on the origin of the inflaton. Modular inflation deals with the inflaton living in 
the closed string sector. On the other hand, brane inflation [14] deals with that in the open 
string sector. In the first case it is sufficient to identify one or more moduli fields which 
are already present in the stabilized compactification scheme and which are displaced from 
the minimum of the potential. In the second case, on the top of the setting required by the 
stabilized compactification, the inflaton field must be engineered by including also probe 
D-branes. The multitude of possible inflationary scenarios in string theory may, at first 
glance, seem confused and far from unique. However, this multitude may all in fact be 
realized if the paradigm of the string theory landscape is adopted |l5| , leading to a picture 
where different types of inflation may proceed in different regions of the landscape. Indeed, 
even simple field theory models may admit a similar landscape of inflationary possibilities. 

String theory inflation models offer a unique opportunity to study (p)reheating in 
an ultra-violet (UV) complete setting where it is conceivable to determine all couplings 
between the inflaton and the standard model (SM) sector from first principles, rather than 
simply assuming some ad hoc couplings on a phenomenological basis. Moreover, because 
post-inflationary dynamics are extremely sensitive to the details of these couplings, it 
follows that any observables generated during preheating (such as nongaussianities [16, 17] 
or gravitational waves |l8| ) have the potential to provide a rare observational window into 
stringy physics. 

Reheating in brane inflation was studied in detail in a number of papers [||, ||, 19. 2C , 21 



|22fl . One of the most interesting realizations of brane inflation is the "warped" KKLMMT 
model pq| , constructed in the context of the KKLT stabilized vacuum 0. The endpoint of 
inflation in this scenario is the annihilation of a brane-antibrane pair. (The inhomogeneous 
dynamics of this annihilation have been discussed in [24, pfifl.) The details of reheating in 



this model are very complicated and involve several stages of energy cascading, first from 
the D — D pair annihilation into closed string modes, next the closed string loops decay to 
excitations of Kaluza-Klein (KK) modes which finally decay to excitations of open string 
modes on the standard model (SM) brane(s). An adequate description of this process 
requires input from string theory. 

Reheating after closed string inflation, on the other hand, is usually assumed to occur 
in the regime where ordinary quantum field theory (QFT) is applicable. In this paper we 
investigate in detail the theory of reheating after closed string modular inflation [26, 27] 
models. We will focus our attention in particular on the scenario of Kahler moduli [26] 
or Roulette |27] inflation models based on the Large Volume Compactification scheme 
of ||, ^, [l0|]. In this model the role of the inflaton is played by a Kahler modulus, t, 
(corresponding to the volume of a 4-cycle of the internal Calabi Yau manifold) and also by 
its axionic partner, 9. However, some of our results may be applicable also in other modular 
inflation models, such as racetrack inflation 2J] based on the KKLT compactification [^]. 

We find that reheating after modular inflation can be quite multifarious and may 
proceed through a variety of different channels (including both perturbative and nonper- 
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turbative effects and also intrinsically stringy physics). The precise identification of these 
decay channels may depend on model building details, such as the location of the SM in 
the Calabi Yau (CY) compactification manifold. In all cases, however, the initial stages 
of the decay of the inflaton in modular inflation proceed via very strong nonperturbative 
preheating decay channels. This is due to the specific shape of the effective inflaton poten- 
tial V(t), which is very nonlinear and which has a sharp minimum. Preheating proceeds 
through a combination of both tachyonic (spinodal) instability and broad-band parametric 
resonance and leads to the copious production of 5t inhomogeneities (particles). Within 
2-3 oscillations of the background field, nonlinear backreaction effects become important 
and the homogeneous condensate is completely destroyed. In fact, this is perhaps the most 
violent known example of preheating after inflation! 

In order to understand the full dynamics of reheating in modular inflation we must also 
identify the decay channels of the inflaton into the visible SM sector. In the case where the 
SM is incorporated on a D7 brane wrapped on the inflationary 4-cycle, the Kahler modulus 
r decays via a direct coupling to brane-bound SM gauge bosons. The initial stages of this 
decay are nonperturbative and involve parametric resonance of the gauge fields while the 
later stages involve perturbative decays. We also consider interactions between the inflaton 
r and brane-bound MSSM fermions such as the Higgsino and gaugino. 

On the other hand, the D7 brane construction described above may result in dangerous 
^-corrections to the inflaton potential which violate the smallness of the slow roll param- 
eters and spoil inflation p9|| . Therefore, it may be desirable to exclude such a wrapping. 
In this case the SM can still be localized on a D7 wrapping some non-inflationary 4-cycle 
of the CY compactification. Such a configuration forbids any direct coupling between the 
inflationary sector and the SM sector, thus complicating the process of reheating. In this 
case the inflaton may still decay to SM states via a nontrivial mixing between the inflaton 
fluctuations and the fluctuations of the moduli associated with the 4-cycle that the SM D7 
wraps. The latter may couple directly to brane-bound SM states. Since the intermediate 
stages of the inflaton decay involve bulk states, the mixing proceeds via Planck suppressed 
operators. 

Finally, we have identified another reheating mechanism which involves distinctly 
stringy physics and which does not require a D7 brane to wrap the inflationary 4-cycle. 
An important model parameter is the value of the inflationary 4-cycles volume, r, at the 
stable minimum of the effective potential, T m . If r m > l s then the supergravity description 



remains valid during both inflation and reheating [26]. In this case, reheating treatment 
is purely field-theoretical and involves calculating the couplings between the inflaton and 
standard model degrees of freedom, as described above. However, as was noted in pTfl , 
although the choice of r m does not alter the field-theoretical treatment of inflation, it does 
crucially impact the dynamics of reheating. 1 If r m j n < l s then the supergravity approxi- 
mation is valid only during the slow roll (large r) regime. At the small values of r relevant 



1 A simple field theoretic analogue of this scenario is a toy model where the inflaton potential is extremely 
flat (to provide sufficient inflation) with an extremely steep minimum. For this toy example ordinary QFT 
is valid during the inflationary stage. However, in the limit that the mass at the minimum approaches m s 
the field theoretical treatment of reheating breaks down and one must instead turn to a stringy description. 
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for (p)reheating stringy degrees of freedom, in addition to the supergravity ones, will be 
excited. In this case the physics of reheating will change drastically from the naive picture. 
One expects, along the lines of ]|, that when r becomes of order the string scale, light 
winding modes are created, and these subsequently decay into free closed strings which cas- 
cade into KK excitations. These intermediate KK modes can, finally, decay into SM states 
on the brane as in ]5|, §]. One can also think about the shrinking 4-cycle as an enhanced 
symmetry point, associated with the quantum production of light degrees of freedom p0[ . 

This paper is organized as follows. In section [2] we review the large volume compactifi- 
cation of type IIB string theory and also discuss Kahler Moduli/Roulette inflation models 
embedded within that setting. In section ||| we discuss the decay of the inflaton in mod- 
ular inflation via self interactions, showing that fluctuations of the Kahler modulus r are 
produced copiously at the end of inflation via strong preheating effects. We study this 
explosive particle production first by solving the linearized equations of motion for the r, 
9 fluctuations and next by performing fully nonlinear lattice field theory simulations. In 
section |] we discuss the decay of the inflaton into SM particles (both perturbative and 
nonperturbative) for the scenario where the SM lives on a D7 wrapping the inflationary 
4-cycle. In section [| we discuss the decay of the inflaton in SM particles for the case 
where the SM D7 instead wraps some non-inflationary 4-cycle of the CY. In section ||] we 
propose another mechanism for reheating in modular inflation which involves distinctly 
stringy excitations. Finally, in section 0, we summarize our results and conclude. 

2. Kahler Moduli/Roulette Inflation in the Large Volume Compactifica- 
tion 

In this section we briefly review recent progress in constructing stabilized "large volume" 
compactifications in type IIB string theory ||, |9j, [[(J and also describe recent efforts to 



embed inflation into such constructions [26, 27] with the role of the inflaton played by the 



Kahler modulus of the internal CY manifold and its axionic partner. 
2.1 Large Volume Compactification 

Let us first discuss the "large- volume" moduli stabilization mechanism of [||, ||, [l0|]. In 
this scenario the Kahler moduli of the CY manifold are stabilized by both perturbative 
and non-perturbative effects. As argued in H, 0, 1C], a minimum of the moduli potential 



in the corresponding effective Ad theory exists for a large class of models. 2 An effective 
\d N = 1 supergravity is completely specified by a Kahler potential, superpotential and 
gauge kinetic function. In the scalar field sector of the theory the action is 



5V=i = Id x 



M 2 

where 



(2.1) 



V(((>, <j>) = e KIM P [IC^DiWDiW - tptWW) + D-terms. (2.2) 



Ml J 



2 The only restrictions being the existence of at least one blow-up mode resolving a point-like singularity 
and also the requirement that ft, 1 ' 2 > ft, ^ > 1, where ft 1,2 , ft 1 ' 1 are the Hodge numbers of the CY. 
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Here /C and W are the Kahler potential and the superpotential respectively, Mp is the 
reduced Planck mass Mp = l/V8irG = 2.4 x 10 18 GeV and (f> 1 represent all scalar moduli. 

The a /3 -corrected Kahler potential [31], after stabilization of the complex structure 
and dilaton, is 

/V 21n(V + |] + ln(<? s )+/C cs , (2.3) 



Ml \ 2, 

Here KL CS is some constant, V is the volume of the CY manifold M in units of the string 
length l s = 27rv / c7 and g s is the string coupling. The second term ^ in the logarithm 
represents the a'-corrections with £ = — ^2(2-1^ proportional to the Euler characteristic 



x(M) of the manifold M. The Kahler metric appearing in (2T) was computed explicitly 
in |7|. 

The superpotential depends explicitly upon the Kahler moduli Tj when non-perturbative 
corrections are included 

W = 9 ^ (w + h ^A ie - a ^ . (2.4) 

Here Wq is the tree level flux-induced superpotential. The exponential terms Aie~ a * T * arise 
due to non-perturbative (instanton) effects such as a gaugino condensate on the world- 
volume of a D7 brane, or Euclidean D3 brane instanton (see e.g. |7|, [||, |9], |l0|). The Kahler 
moduli are complex, 

Ti = n + iOi , (2.5) 
with Ti the volume of the z-th 4-cycle and 6i its associated axionic partner, arising from 



the Ramond-Ramond 4- form C4. In (2^) the constants A^, ai depend upon the specific 
nature of the dominant non-perturbative mechanism. For example, Oj = y- for Euclidean 
D3-brane instantons. 

In the simplest cases (such as 1 1 g g j) the volume V can be written in terms of the 
Ti as follows: 

V = a^-±X l rf^ . (2.6) 

Here a and Aj are positive constants depending on the particular model. For example, 
the two-Kahler model with the orientifold of 1169 ] studied in || 10, |33j has n = 2, 
a = l/9\/2 and A 2 = 1. 



The formula ( |2.6| ) suggests a "Swiss-cheese" picture of a CY, in which t\ describes 
the 4-cycle of maximal size and controls the overall size of the CY. This modulus may be 
arbitrarily large. On the other hand, T2, . . . ,r n , are the blow-up cycles which control the 
size of the holes in the CY. These moduli cannot be larger than the overall size of the 
compactification manifold. This CY manifold is schematically illustrated in Fig. |]. 

Including both leading perturbative and non-perturbative corrections one obtains a 
potential for the Kahler moduli which, in general, has two types of minima. The first type 
is the KKLT minima . Relevant for the present study are the "large- volume" AdS minima 
studied in [|, |, 0. These minima exist in a broad class of models and at arbitrary values 
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Figure 1: A cartoon of the large volume compactification manifold illustrating the ingredients 
relevant for Kahler moduli inflation and (p)reheating. The modulus t\ = Re(Ti) controls the 
overall size of the compactification while the moduli Tj = Re(Ti) [i > 2) control the size of the 
blow-up 4-cycles (hole sizes). We have labeled as the last 4-cycle to stabilize, and hence this 
modulus is associated with the final observable phase of inflation. The total volume t\ and other 
cycles (for example T3) are assumed to already be stabilized. In the text we consider two possible 
scenarios for the location of the SM: (i) a D7-brane wrapping the inflationary 4-cycle, T2, and, 
(ii) a D7-brane wrapping the non- inflationary 4-cycle, T3. We have illustrated these wrappings 
schematically. 



of parameters. An important characteristic feature of these models is that the stabilized 
volume of the internal manifold is exponentially large, V m ~ exp (ar m ) (here V m , r m denote 
the values of V and T2 at the minimum of the potential), and can beO(10 5 - 10 20 ) in string 
units. The relation between the Planck scale and string scale is 



(2.7) 



Thus, these models can have m s in the range between the GUT and TeV scale. 



As argued in ||, 10], for generic values of the parameters Wq, A^, ai one finds that 
T i 3* Tj (i > 2) and V S> 1 at the minimum of the effective potential. In other words, the 
sizes of the holes are generically much smaller than the overall size of the CY. 

An important consistency condition for the supergravity approximation to be valid is 
that the value of each Tj at the minimum of the potential should be larger than a few. 
This criterion ensures that the geometrical sizes of the 4-cycles of the CY are all larger 
than the string scale. One expects any violation of this condition to be associated with the 
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production of stringy (as opposed to field theoretic) degrees of freedom, such as winding 
modes. More on this later. 



2.2 Roulette Inflation 

Let us now consider inflation in the context of the large volume compactification described 



above 27]. The scenario we have in mind is the following. Suppose all the moduli Tj 
are initially displaced from their minima. The dynamics will then drive the various fields 
Tj to roll towards their respective minima. For inflation, we focus on the last modulus 
T 2 = t + i9 to reach its minimum so that r, 9 are still dynamical while all other moduli 
(in particular the total volume V) are stabilized. 3 See Fig. |l] for a cartoon of this scenario. 
Note that in the Roulette inflation model both fields r and 9 play an important role in 
driving inflation. After having fixed all moduli Tj (i ^ 2) we find an effective potential for 
the inflaton fields r, 9 of the form 

(2.8) 



V(r,9) = g 



g 4 s \8(a 2 A 2 ) 2 ^e- 2a2T 4W a 2 A 2 Te- a * T cos (a 2 9) Ay 



3aX 2 V m V 2 



where we have expanded to order in the (exponentially) small parameter 1/V m . In 

J2.8D the uplifting Ay is some model-dependent constant which should be tuned so that 
V = at the minimum. 

The potential surface V(t, 9) has a rich structure, illustrated in Fig. 0. The potential is 
periodic in the axion direction and is exponentially flat (but slightly rippled) at large values 
of the radial coordinate r in field space. The minimum (near r = 0(1) and cos(a2#) = — 1) 
is extremely steep. At r —* there is a sharp potential barrier, which we have cut-off 
in Fig. ^ in order to make the salient features of the potential more transparent. The 
shape of the potential near the edge of the minimum is locally reminiscent of the the 
racetrack inflation potential |28[| , as well as the "natural inflation" potential involving a 
pseudo-Goldstone boson |}2| (except that in our case both 9 and r must be simultaneously 
considered). 

The potential ( |2.£| ) allows for a rich ensemble of inflationary trajectories in the Kahler 
moduli space p7|] , depending on the choice of initial conditions and model parameters. 
We have superposed several of these trajectories on the potential surface in Fig. [2] for 
illustration. Inflationary trajectories may undergo many revolutions at roughly constant 
r, spiraling along the angular direction, 9, like a ball on a roulette wheel. At some point, 
the inflaton eventually gets caught in the trough along 9 = (where I is an integer) 

and rolls down to small r. Generically the last few e-foldings of inflation occur along this 
trough, which is stable against 59 perturbations and along which V(r, 9) takes the form 
considered in [p^j . Inflation ends when the inflaton rolls into the sharp minimum, like the 
ball on a roulette wheel falling into the pocket. 

Models of this type necessarily introduce a statistical ("gambling") aspect to the con- 



straints that observational data imposes on inflationary model building. In 27] it was 



3 The consistency of this approach was discussed in E^], see |m| for more general types of inflationary 
trajectories. 
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Figure 2: The T2 potential V(t, 9) surface for a representative choice of parameters, using polar 
coordinates to illustrate the periodic structure of the potential in the axion direction, 6. Super- 
imposed on the potential surface are three different inflationary trajectories, showing the rich set 
of possibilities in Roulette inflation. Inflation proceeds in the large r region where the potential is 
exponentially flat. On the other hand, preheating after inflation takes place during the phase of 
oscillations the extremely steep minimum near r = 0(1) and cos (02 9) = — 1. 



advocated to view this theory prior as a probability distribution on an energy landscape 
with a huge number of potential minima. In |3£| the spectrum and nongaussianity of 
the curvature fluctuations for the variety of Roulette inflation trajectories was discussed. 
In [39] cosmological fluctuations were studied more generally in multi- field models with 
non-standard kinetic terms and the Roulette model was considered as a special case. 

As discussed previously, consistency of the supergravity approximation requires ad- 
justing the parameters so that r m > a few. We note, however, that even if the SUGRA 
description in terms of the scalar potential is not valid at the minimum, it still can be 
valid at large r, exactly where we wish to realize inflation. Hence a small value of T^n 
does not constrain the Kahler Moduli/Roulette inflation. Rather, the only consequence of 
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small T m is that the endpoint of inflation, reheating, would have to be described by string 
theory degrees of freedom (rather than simple field theoretic objects). In this paper, we 
will propose one of the possible scenarios of this type. 

Finally, we note that there are several types of perturbative corrections that could 
modify the classical potential (p.S| ) on Kahler moduli space: those related to higher string 
modes, or a'-corrections (coming from the higher derivative terms in both bulk and brane 
effective actions) and also string loop, or ^-corrections (coming from closed and open string 
loop diagrams). The most dangerous corrections that could spoil exponential flatness of the 
potential at large t, are the latter category: those coming from open string diagrams. This 
type of corrections are relevant for the models where the SM lives on a D7 brane wrapped 
on the 4-cycle associated with the inflationary modulus, 4 t 2 , and are expected to spoil the 
flatness of the inflaton potential; see [^] for an estimate of this effect. One can evade this 
problem simply by excluding such a wrapping and assuming that the SM lives on a D7 
wrapping some other (non-inflationary) hole in the compactification. 5 This scenario has 
the disadvantage of complicating the reheating process, since it forbids a direct coupling 
between the inflaton and the visible sector. In this paper, we proceed phenomenologically 
and consider models both with and without a D7 brane wrapped on the inflationary cycle. 
In both cases we identify the dominant decay channels of the inflaton into SM particles. 

It will be handy for the discussion below to have estimates of the masses of the fields 
playing a role in inflation and reheating. Recall that from equation fl2.7|) the string scale 
and Planck scale are related as m s /M p ~ V -1 / 2 . The gravitino mass is m 3/ / 2 /Mp ~ V -1 
and the masses associated with the deformations of the total volume and hole size(s) are, 
respectively, given by m Tl /M p ~ V _3/2 (ln V)' 1/2 , m n /M p ~ InV/V (for i > 2). See §7| 
for more detailed discussion. 



3. Preheating via Self-interactions in the Inflaton Sector 



In this section we study the energy transfer from the homogeneous inflaton condensate into 
fluctuations. We proceed by first considering the dynamics of linear perturbations about 
the homogeneous inflaton background and next by studying the fully nonlinear dynamics 
of the system using lattice field theory simulations. In the next section we study the 
subsequent transfer of energy from these fluctuations into excitations of the SM particles. 

3.1 Equations for Linear Fluctuations 

We consider the endpoint of inflation in the Roulette inflation model reviewed in the last 
section. To this end we study the action 



S 



-9 



l -n - * /c 22 (d^r + d^e) - v(r, o) 



2 2 

where the (2, 2) component of the Kahler metric depends only on the 4-cycle volume 

3aA 2 [2V m + f + 6aA 2 r 3 / 2 ] _ 3aA 2 1 



/C 



22 



2(2V m + 2 V^ 



4V m y/r 



+ • 



(3.1) 



(3.2) 



4 This possibility is illustrated schematically as "scenario 1" in Fig. 
5 This possibility is illustrated schematically as "scenario 2" in Fig. 
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The potential V(r, 9) is given explicitly by (2^8) and in the second equality of fl3.2|) we have 
expanded to leading order in V^ 1 . Throughout the rest of this paper we use Planck units, 
Mp = 1, although we occasionally write out the factors of M p explicitly. For simplicity 
we restrict ourselves to "parameter set 1" (as defined in |27j) corresponding to the choice 
Wo = 300, a 2 = 2vr/3, A 2 = 0.1, A 2 = 1, a = l/(9\/2), £ = 0.5, g s = 0.1 and V = 10 6 . We 
have considered also different parameter choices and found the results to be qualitatively 
similar in all cases. 

To understand preheating in the model ( |3.1[) , let us first consider the linearized equa- 
tions for the fluctuations. In order to make the analysis tractable we introduce the canonical 
inflaton field 



defined by dtp = ±y/fC 22 dr so that <fi = W f^p^ 2 - t 3//4 in the large volume 
limit and with some arbitrary convention for the origin of field space. Now the action (^ 
takes the form 



S 



I 



d xyj—g 



\{d<j>f 



1 



+ -<?M{doy + u{<j> i o) 



where 



U{4>,6) 



-ln[/C 25 (r(0))] 
V[r(ct>),6] 



(3.3) 

(3.4) 
(3.5) 



Our choice of notation follows |39| . 

For generic inflationary trajectories, the final few e- foldings of inflation take place 
along the trough 9 = ir(2l + l)/a 2 (with I integer) |2"^| . Thus, to study the linear regime 
of particle production during preheating after inflation we expand the fields as 



0(t,x) 
9(t,x) 



0o(t) + ^(t,x) 
(21 + l)vr 



09 



(3.6) 
(3.7) 



In the left panel of Fig. ^ we plot the effective potential for the homogeneous motion of the 
canonical inflaton <f> along the 9 = (21 + l)ir/a 2 trough, that is U [(/>, (21 + l)ir/a 2 ]. This 
potential displays a long flat region at large eft which is relevant for inflation and a very 
steep minimum at (f> = 4> m . Preheating during inflation takes place during the phase of 
oscillations about this minimum. In the right panel of Fig. |3| we plot the time dependence 
of the inflaton condensate (j>o(t) during this oscillatory phase. Due to the extreme sharpness 
of the potential minimum the inflaton passes very quickly through the region close to <p m 
while the extreme flatness at larger <f> means that the inflaton spends a long time near the 
right-hand side of the valley. We denote the period of oscillations about the minimum by 
T and note that for typical parameters T ~ 10~ S H~ 1 . 

Let us now consider the linear fluctuations 6(f), 69 in the inflaton and axion fields. The 
equations of motion for the Fourier modes 5(pk and 59k take the form 

Ik 2 

b k + ZH8(f> k + 



S9 k + [3H + 2b^} 59 k 



k 2 



+ e~ 2b U 



59 k 











(3.8) 
(3.9) 
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Figure 3: The left panel shows the effective potential for the canonical Kahler modulus 4> along the 
axion trough, U((f>, (21 + 1)77/02) showing the long exponentially flat region relevant for inflation 
and also the steep minimum relevant for the preheating phase of post-inflationary oscillations. We 
have labeled the point where inflation ends (where the e slow roll parameter is unity) and also the 
point where the effective mass-squared V"((j>) flips sign, corresponding to the cross-over between 
the tachyonic and non-tachyonic regions. The right panel shows the oscillatory time evolution of 
the homogeneous inflaton <f>o(t) at the end of inflation. 



where the quantities U(4>, 9) and b(4>) (and their derivatives) are understood to be evaluated 
on the unperturbed values <fi = (f>o(t) and 9 = 9o(t) = (2Z+1)tt / 'a%. Neglecting the expansion 
of the universe 6 and defining the canonical axion fluctuation 8ij) as 



5^ k = e b 59 k 

we obtain oscillator-like equations for the mode functions 



5ii> k + vl k 8ipk 








In ( 3.11 , 3.12 ) the effective frequencies are 



2 



2 



k + U M 
k 2 + Ml cS (t) 
k 2 - b t< p(f> - (i 
k 2 + Ml cS (t) 



+ b 



(3.10) 



(3.11) 
(3.12) 



(3.13) 



(3.14) 



and we have introduced the notation M 2 S (t) for the time-dependent effective masses of the 
fields. We solve equations ( 3.11 , 3.12 ) numerically. As usual Q we initialize the modes as 
54>k = l/y^u^.fc, 5cpk = — ^v^^fe/2 at t = (and similarly for 5ipk) which, physically, 
corresponds to starting with pure quantum vacuum fluctuations with occupation number 



°This neglect is justified since typically the total time for preheating to complete is of order 10 3 H 1 . 
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rik = 0. Violations of the adiabaticity condition |cjfe|/a;| <C 1 are associated with particle 
production and lead to the generation of nonzero occupation number, n k 7^ 0. As is typical 
in the theory of preheating we define the occupation number for the inflaton fluctuations 

as 

_ u<j>,k 



+ \S<f> k \ 



\ (3-15) 



and similarly for the occupation number nt associated with the axion. 
3.2 Instability of Kahler Modulus Fluctuations 

Let us first consider the equation (3.11) for the fluctuations 54>k of the (canonical) inflaton. 



This equation can be viewed as an effective Schrodinger-like oscillator equation with "po- 
tential" determined by the effective mass term, M 2 eS (t). In Fig. || we plot the behaviour 
of this effective mass (the dashed blue curve) as a function of time. We see that M? eff has 
a very particular time dependence: the fluctuations 5<p are nearly massless during the infla- 
ton oscillations except for periodic "spikey" features uniformly spaced at intervals At = T. 
As long as the adiabatic invariant lu^l/w^ remains small, |wfe|/a;? <C 1, particles are not 
produced and the occupation number n k will be close to a constant. For modes with wave- 
number k 2 sufficiently small, the spikey structure in M 2 eS leads to extreme violations of 
adiabaticity: |cjfc|/u;| ^> 1. In Fig. ^ we plot the adiabaticity parameter taking k = for 
illustration (the dashed blue curve). Notice that each spike in M 2 eS is also accompanied 
by a tachyonic phase where M? ff < 0. 

The periodic time-dependent behaviour of the effective mass M? eff leads to unstable 
momentum bands where the modes 5(f)f~ grow exponentially as 

6<f> k (t) ~ e^ T f k (t/T) (3.16) 

with fk{t/T) some periodic oscillatory function. We computed the Floquet exponent p, k 
appearing in (3.16) numerically and the result is displayed in Fig. 0. The broad unstable 



region at low k corresponds to tachyonic growth of IR modes coming from the Mj eff < 
regions in Fig. ^[ On the other hand, the UV region of Fig. [6| featuring a band structure is a 
result of broad-band parametric resonance. The modes S(f) k belonging to these two regions 
of phase space display very different behaviour and both tachyonic and resonance effects 
play a crucial role in the dynamics of preheating. In Fig. ^ we plot the time evolution of 
the Kahler modulus mode functions for k = 0.01 (corresponding to the IR tachyonic 
regime), k = 0.5 (corresponding to the UV regime of broad-band parametric resonance) 
and k = 0.08 (corresponding to the intermediate regime where these two effects cannot 
be disentangled). Here denotes the mass of the canonical inflaton at the minimum of 
the potential. Note that these different kinds of preheating can only be separated in the 
linear theory; at the nonperturbative level such a distinction is impossible. 

To illustrate how violent the process of preheating in modular inflation is, we consider 
the energy density in a given wave-number k, which coincides with ~ /c 4 n^ if we evaluate 
Pk at the point in the inflation oscillations where Ma, e g = (this point is illustrated on 
the potential in the left panel of Fig. 0). This quantity is plotted as a function of k in 
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Figure 4: The effective mass-squared of both the Kahler modulus <j> (the dashed blue curve) and 
axion ip (the solid red curve) as a function of time. We plot the effective mass in units of m<t,, the 
mass at the minimum of the potential. Note that M e g(t) actually exceeds during the inflaton 
oscillations. This corresponds to the steep curvature on the left-hand-side of the potential minimum 
in Fig. §. 



Fig. |8] for three different time steps in the (linear) evolution. We see that within only three 
oscillations of the background field (that is, at t = 3T) the fluctuations contain many orders 
of magnitude more energy than the condensate. At this point backreaction becomes critical 
and the linearized treatment breaks down. Therefore, in order to study the dynamics after 
the first 2-3 oscillations one must turn to nonlinear lattice simulations, which we consider 
in the next subsection. 

We have also studied the dynamics of the linear fluctuations of the axion, equation 
p. 12 ). The effective mass M| eff (i) for the fluctuations 5ipk of the (canonical) axion. 
is plotted in Fig. |4| (the solid red curve). The behaviour is qualitatively similar to the 
inflaton effective mass, however, the spikes are less sharp and the tachyonic regions are 
absent. Thus, we expect the production of axion fluctuations to be less efficient. This 
intuition is confirmed in Fig. [5] where we show that adiabaticity is never violated for these 
modes (see the solid red curve). Therefore the endpoint of modular inflation is dominated 
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Figure 5: The quantity \dM c s/dt\/\M^ s \ for both the Kahler modulus (j> (the dashed blue curve) 
and axion ip (the solid red curve) as a function of time. This quantity provides a measure of the 
violation of adiabaticity and coincides with the adiabatic invariant ujk/ojf, in the IR. The spikey 
structure of the Kahler modulus effective mass leads to extremely strong violations of adiabaticity 
and particle production. On the other hand, the axion effective frequency varies slowly during the 
inflaton oscillations and axion particles are not produced. 

by the extremely nonperturbative production of Kahler modulus fluctuations, rather than 
their axionic partner. During later stages of reheating, however, axions may be produced 
by rescattering effects. 

3.3 Lattice Simulations of Preheating in Roulette Inflation 

In the last subsection we studied the dynamics of the linear fluctuations St, 59 at the 
endpoint of inflation. We found that the particular shape of the inflaton potential (which 
is highly nonlinear and has a very steep minimum) imparts an extremely non-adiabatic 
time dependence to the 4-cycle modulus r and leads to explosive production of 5t particles. 
Hence the linearized analysis of this dynamics breaks down very rapidly, within just 2-3 
oscillations of the background field, and a full treatment requires fully nonlinear lattice 
field theory simulations. Hence we now study the evolution of the fields r, 8 in a fully 
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Figure 6: The dimensionless characteristic exponent Hk (Floquet exponent) for the exponetial 
instability of the (canonical) Kahler modulus fluctuations S<f>k, defined in ( 3.16 ). The broad unstable 
region in the IR comes from the tachyonic regions where Mj eff < whereas the UV region displays 
the band structure that is characteristic of parametric resonance. The behaviour of the modes 5<pk 
in these two regions is qualitatively different. 



nonlinear way, taking into account also the expansion of the universe self-consistently. 

Our numerical simulations are done using a new lattice code which contains (modified) 
elements of both LatticeEasy [35] and DEFROST The first version of this code was 



employed in the paper [ 40 1 to study infra-red cascading during inflation. Here we have 
further generalized the code to allow for scalar fields with non-canonical kinetic terms. 
With the MPI-parallelized code we have used 64 processors on the CITA cluster to evolve 
the fields r and 9 in a cubic box with 512 3 grids. Vacuum mode functions are put in as 
initial conditions. The very first stage of the evolution, when inhomogeneity can be treated 
linearly 7 , is performed in momentum space. From this initial stage of momentum-space 
evolution we can determine which Fourier modes r^, 9^ experience unstable growth during 
preheating and ensure that the relevant dynamical scales are captured by our simulations 
(that is, the /c-modes which dominate the unstable growth are well between the IR cut-off 



7 In the case of Roulette inflation this occurs only during inflation and the first oscillation of preheating. 
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Figure 7: The behaviour of linear Kahler fluctuations during preheating after modular inflation 
illustrating the combination of tachyonic instability and parametric resonance. The left panel 
shows the mode behaviour for k — O.Olm^ (corresponding to the IR tachyonic regime) while the 
right panel shows the mode behaviour for k = 0.57710 (corresponding to the UV regime of parametric 
resonance). The middle panel is k = O.O8m0, corresponding to the intermediate regime where both 
effects are active. 

k[R = L" 1 denned by the box size L and the UV cut-off kuv = Ax -1 denned by the 
lattice spacing Ax). Once the particle occupation numbers become large we switch over to 
configuration-space evolution (at the switching time the gradient energy is roughly 10 
as compared to the energy in the homogeneous condensate). We subsequently run our 
configuration-space evolution for a total time At = 0(3 — 4)T (where T is the period 
of oscillations of the homogeneous background). This time scale is more than sufficient 
to see the homogeneous inflaton condensate completely decohered into inhomogeneous 
fluctuations. 

To illustrate the rapid and violent development of inhomogeneities during modular 
preheating we study the energy density in gradients of the Kahler modulus r, given by 

Pgrad = JC jTT ^diTd l T = J ^ K-,tt \ &T k \ 2 (3.17) 

In Fig. |{| we plot the power in gradient energy, given by /c 5 /C iTT |<5rfc| 2 /(27r 2 a 2 ), for three time 
steps in the evolution. After the first oscillation, at t ~ T, we see the band structure of 
growing modes in the UV is characteristic of parametric resonance. However, by the time 
t ~ 2T the evolution has become completely nonlinear leading to the destruction of the 
band structure and the cascading of power into both the IR and UV. By the time t ~ 3T the 
energy density in gradients is comparable to the homogeneous background energy density, 
illustrating the decay of the condensate. By this time, however, our numerical evolution 
has become unreliable since the power in inhomogeneous fluctuations has cascaded down 
to the UV cut-off, /ctjv = Ax" 1 ; beyond this point we are bleeding energy into the grid. In 
any case, at later times decays into fields which were not included, such as brane-bound 
gauge fields and fermions, start to become important. We consider these decays for a 
variety of scenarios in the remainder of this paper. 

We have also considered the dynamics of the expansion of the universe during pre- 
heating. The effective equation of state during preheating is to = P/p = 0.1, which is very 
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Figure 8: The energy spectrum ~ k A nk for St, calculated with linear theory, showing the rapid and 
violent production of particles after modular inflation. Within three oscillations of the background 
field (by t — 3T) the energy density is significantly larger than the energy in the homogeneous 
condensate and at this point nonlinear feedback must be taken into account. 

close to a matter dominated expansion. This is consistent with Fig. || where we see that, 
in the linear theory, the 5t particles which dominate the energy density of the universe 
have momentum in the band kjm^ ~ 0(0.01 — 0.1) and hence we expect these to be non- 
relativistic near the minimum: u>k = ^jk 2 + = m$ + + • • • . (Note also that the 
Floquet exponent jjL^ in Fig. || peaks over roughly the same range of scales.) 

The violent production of inhomogeneities leads to significant backreaction effects 
which alter the form of the effective potential, in particular the location of the minimum 
and the effective mass at that minimum, m T . The late-time perturbative decays of the 
produced inflaton fluctuations 5t into SM fields (which we study in detail in the remainder 
of this paper) depend sensitively on the value of ra T at the minimum, hence it is impor- 
tant to determine whether this differs significantly from the value one obtains neglecting 
backreaction effects. 

To illustrate how backreaction effects modify the effective potential, consider a simple 
scalar field with equation of motion n<j> = V (</>). Writing (j) = ((f)) + 5(f) and taking the 
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Figure 9: The power in Kahlcr modulus gradient energy, which is given by /c 5 /C jrr |5Tfc| 2 /(27r 2 a 2 ), 
from our lattice simulations, normalized this to the total background energy density, p. The spec- 
trum is shown for four time steps in the evolution: t = T,t = 2T, t = 2.25T and t = 3T, illustrating 
the rapid decay of the homogeneous inflaton condensate into inhomogeneous fluctuations. See the 
text for further discussion. 

average of the Klein-Gordon equation we find an equation for the mean field 

□<^> = (V l(4>) + 64>}) 

= [(</>)] (3.18) 

In general the effective potential V e s will differ from V(<j>). This can be illustrated explicitly 

<s</> 2 

by defining the averaging procedure as (V ((/))) = —r= — J d[8(f)]e 2,T <t>V' '{(f)) with the 
variance of 4>. In general this averaged force differs from V' (</>). Note that in our case both 
inhomogeneous fluctuations and nonlinear effects are large, thus a Gaussian distribution 
for 5t = t — (t) may not be appropriate. However, this does not alter the qualitative 
argument that we are making. 

We have estimated the effective mass for the r fluctuations at the minimum, m T , in- 
cluding backreaction effects, by extracting the rate of oscillation of the (spatially) averaged 
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field (r) from our lattice simulations. We find that 

m^ = j3mf d (3.19) 

with j3 ~ 2.5 for our choice of parameters. Hence, the true mass including backreaction 
effects differs from the mass without backreaction by a factor order unity. However, since 
the decay of St into brane-bound gauge fields depends on the effective mass as T oc rr? T the 
value of the factor (3 may be important. 

We have also investigated the production of gravitational waves from the rapid devel- 
opment of inhomogeneities associated with the violent decay of the inflaton after modular 
inflation. We have quantified the production of gravitational waves numerically, following 
the approach adopted in [18]. We did not find any significant production of gravitational 
waves for the model parameters considered. 

4. Transfer of Energy into the SM Sector: D7 Wrapping the Inflationary 
4-Cycle 

We have seen that the initial stages of preheating in Roulette inflation proceed nonper- 
turbatively by a very strong instability and lead to explosive production of St fluctuations 
(particles) which destroy the homogeneous inflaton within just a few oscillations. How- 
ever, this violent particle production is not sufficient to make contact with the usual hot 
big bang picture. We must also ensure that the energy in St fluctuations can be efficiently 
transferred into excitations of the visible (SM) sector. Determining the details of this en- 
ergy transfer requires a complete understanding of how the inflaton fields r, 9 couple to 
SM degrees of freedom and hence this is is necessarily a model dependent issue. Below we 
will proceed phenomenologically and consider a variety of possible scenarios, identifying 
the dominant decay channel of the inflaton in each case. Throughout this section (and the 
next) we assume that the value of r at the minimum, r m , is sufficiently large to validate 
the effective supergravity treatment and hence the details of reheating can be studied using 
only QFT methods. In a subsequent section we will relax this assumption. 

First let us consider the case where the SM lives on a D7 brane wrapping the inflation- 
ary 4-cycle, T2 = r + i6. This wrapping is illustrated schematically in Fig. [I] as "scenario 
1" . From the computational point of view this is the simplest imaginable scenario since it 
allows for a direct coupling between the inflationary r and the fields of the SM. However, 
as discussed previously, such couplings lead to loop corrections to the inflaton potential 
which are expected to spoil slow roll |29|]. Until a robust calculation of such loop effects is 
available we will proceed phenomenologically and suppose that the offending contributions 
to the slow roll parameters can be canceled by some fine tuning (or otherwise). 

4.1 Inflaton Coupling to Photons 

Since the inflaton, r, is the volume of the 4-cycle wrapped by the SM D7 we expect that r 
couples to visible states as an overall prefactor to the SM Lagrangian. Thus, an effective 
Lagrangian for the inflaton-photon interaction takes the form [37] 



£ 7 = ~- a tF^F, v (4.1) 
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where F^ v = d^A u — d v A^ is the field strength associated with the gauge field and 
A = r" 1 once we have noi 
stabilized at the minimum. 



A = r m x once we have normalized A^ so that it has canonical kinetic term when r is 



4.2 Photon Preheating 

The strong nonperturbative production of 5r fluctuations observed in section || provides 
a motivation to look for preheating also in other bosonic fields, such as the photon, A 1 *. 
Therefore we study the stability of the quantum fluctuations of A^ in the background of the 
oscillating homogeneous inflaton To(t). To this end we choose the transverse gauge: ^4° = 0, 
d%A % = 0. The equation of motion for the Fourier transform of the spatial components of 
the gauge field A l (jb, x) in the homogeneous inflaton background takes the form 



A k + 



k 2 

A k + —A k = (4.2) 



a- 



where we suppress the vector index i on the modes Al (t) for ease of presentation. Equation 
( |4.2j ) can be put into oscillator-like form by introducing the field A k = (aro) 1//2 A k . We 
find 

A k = (4.3) 



d 2 2 



where the time- varying frequency is 



2 k 2 ^ \ ( foV lr 1 rr f 1 tt2 1 



% + Ml cE {t) (4.4) 



The time dependence of the effective photon mass M 2 eff (i) (neglecting the expansion of 
the universe) is plotted in the left panel of Fig. [l0|. The behaviour is qualitatively similar 
to the effective mass for the 5t fluctuations (up to an overall sign flip). The oscillatory 
behaviour of the photon effective mass leads to parametric resonance of the photon modes 
A k with stability /instability bands in which the photon fluctuations grow exponentially as 
A k (t) ~ e Mfc '/ T . We have computed the Floquet exponent, [i k) numerically and this result 
is plotted in the right panel of Fig. [R] where we see the characteristic features of narrow 
band parametric resonance. 

We have seen that very early stages of oscillations of the inflaton condensate tq{€) 
leads to copious production of photons living on the SM D7-brane via parametric resonance 
preheating. This phase of violent, nonperturbative photon production is very short-lived 
since the production of 5t inhomogeneities completely destroys the condensate within 2- 
3 oscillations. The later stages of the production of SM photons will instead involving 
perturbative decays r — > 77, to which we now turn our attention. 

4.3 Perturbative Decays to Photons 

After the phase of violent nonperturbative particle production discussed above the homo- 
geneous inflaton r will settle down to the minimum of the potential and reheating will 
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Figure 10: The left panel shows the time dependence of the effective photon mass, M? eff (t). The 



right panel shows the Floquet exponent Hk for the mode functions Ak(t) ~ e 



be dominated by the perturbative decays of produced inflaton fluctuations into SM pho- 
tons. To study such processes let us consider the Lagrangian for the total volume modulus 
n = (where the subscript b stands for "big") and the hole-size ti = t = t s (where the 
subscript s stands for "small"). Writing T{ = (rj) + 5ti (with i = b,s) we have, near the 
vicinity of the minimum 

C = -IC^iSr^iSrj) -V - (M 2 )^(«5r i )(<5T i ) - ^ r s F^F^ + ■■■ (4.5) 

The mass matrix (/C -1 M 2 )ij is not diagonal, nor are the fields 5ri canonically normalized. 
Following [37] we can put this action into a more conventional form by introducing the 
physical modulus fields x defined by 





(4.6) 



\ R s <s> R sx 

The elements of the rotation matrix Rij are computed explicitly in [37], however, paramet- 
rically we have R b $ ~ V 1/6 , R bx ~ V 2/3 , i? s $ ~ V 1/2 and R sx ~ V°. Thus the canonical 
field $ is mostly the small 4-cycle 5t s , while the field x is mostly the big 4-cycle <5r&. In 
terms of the canonical basis x> $ the Lagrangian (4.5) takes the form 

c = -\m 2 - \(dxf -v - ^ - - \f"» Fiw 

+ X X 4>$X$ 2 + ^xx®X 2 + A xxx x 3 + A$*$$ 3 + • • • (4.7) 
The masses of the canonical moduli are 

m - ~ V 3/ 2 (lny)V 2 ^ ^ 
InV 

m$ ~ — M p (4.9) 
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so that <E> is significantly more massive than Xi consistent with the intuition that since 
T b 3> t s the volume modulus should be lighter than the modulus associated with the hole 
size (at the minimum). The light x particles have Planck suppressed couplings to brane- 
bound gauge fields 

X77 MplnV V ' 

since the fluctuations of the total volume are bulk modes. On the other hand, the coupling 
of massive <J> particles to gauge fields are only string suppressed 

The violent production of 5t s fluctuations at the end of inflation leads to a huge number 
of canonical particles near the minimum, in addition to a small admixture of x which 
arise due to the nontrivial mixing of St s and Stj,, see equation ( |4.6[ ). The massive <3? particles 
subsequently decay to photons via perturbative processes such as <I> — > 77 mediated by the 



couplings in (4/7). The decay width was computed in [37] 



*~* 77 ~ 64vr ( } 



Using the estimates (fOj) and (4.11) we find 

(InV) 3 
1 $^ 77 



M p (4.13) 



Also notice that, due to the 8i~b — St s mixing, a certain amount of x particles will be 
produced from the oscillating inflaton field. This does not lead to a cosmological moduli 
problem because the x particles are extremely massive for the value of V which we consider. 
For example, taking V ~ 10 6 — 10 9 we have m x ~ 10 4 — 10 8 GeV which is sufficiently massive 
that the residual x moduli will decay well before Big Bang Nucleosynthesis (BBN) [41]. 8 



4.4 Inflaton Coupling to Fermions 

The inflaton t s will couple not only to photons but also to fermion fields living on the world- 
volume of the D7 brane wrapping the inflationary 4-cycle. The effective action describing 



these interactions generically takes the form [44]: 



£int = (4-14) 

where is the dimensionless coupling and ^> schematically denotes any MSSM fermion. 

The decay rate for this type of interaction is 

h 2 

r^xpvp = — m$ (4.15) 

07T 



8 For very large values of the compactification volume, V ~ 10 9 , the x modulus may decay very close to 
the onset on BBN. For such parameters a more careful treatment of the decay of \> taking into account 
factors of order unity, may be necessary. 
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(in the limit m$ 3> m^). 

Prior to electro-weak symmetry breaking a direct coupling to SM fields (such as the 
electron) is absent and the decay of $ into fermions is dominated by the production of 
Higgsinos — > HH and gauginos $ — ► AA. The relevant couplings were computed in [44] 



h ™ ~ W^v (416) 
*™ ~ W^v (417) 



We therefore have the following estimates for the decay rates 



r *">™ ~ VMnV Mp (418) 
r *-AA ~ Tm^TT^p ( 4 - 19 ) 



We see that the decays of 3> into Higgsinos dominate over the decays into gauginos. 
4.5 Reheating Temperature 

Let us now estimate the reheating temperature for the scenario where the SM lives on 
a D7 wrapping the inflationary cycle T%, Note tlicit I"^— ^--y-y/r^ ^.jjjj 

(InV) 4 ~ 40 4 

for a compactification volume of order V ~ 10 6 — 10 9 , which is of interest for Roulette 
inflation. Thus, comparing equations (4.43) and (4.1 g|) we see that the lifetime of the 



inflaton is dominated by its decay into photons. Using the standard result from the theory 
of reheating || the reheat temperature is given by 



(InV) 



3/2 



T r ~ O.lyTMp ~ O.P — ^ — M p (4.20) 

Using V ~ 40 6 we obtain T r ~ 10 13 GeV and using V ~ 10 9 we have T r ~ 10 10 GeV. 
Normally such a high reheat temperature would lead to a gravitino problem (specifically 
the late decays of the gravitino interfere with the successful predictions of Big Bang Nu- 



cleosynthesis). However, in our case the gravitino mass is extremely high [37|: 



m 3/2 ~ (4.24) 

For V ~ 40 6 — 40 9 we have 7713/2 ~ 40 9 — 40 12 GeV, hence the gravitino is so massive 
that it decays well before nucleosynthesis |4i| . The disadvantage of such a high gravitino 
mass temperature is that SUSY is broken at too high a scale to explain the electro-weak 
hierarchy of the SM. An alternative scenario to explain the fine-tuning of the Higgs mass 
is through vacuum selection effects on the string landscape ]4^| . On the other hand, see 
[43|| for a scenario where this gravitino mass is compatible with TeV scale soft terms. 

There is, however, another cause for concern. At sufficiently high temperatures thermal 
corrections to the effective potential tend to destabilize the moduli and drive them to 
infinity, leading to decompactification || . The decompactification temperature, T max , 
can be accurately approximated by T max = V^ 4 where = m^ 2 M p is the height of 
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the potential barrier separating the large volume AdS minimum from the supersymmetric 
minimum at infinity in moduli space. The maximum allowed reheat temperature (above 
which the internal space decompactifies) is therefore 

rr rv, M P (4.22) 



V 3 / 4 



Comparing (4.2C) to (|4.22j) one sees that, keeping track only of factors of the total 
volume V, we have T r is marginally below the decompactification temperature. However, 
this conclusion could easily be altered by factors of order unity which we have not taken 
into account. Such factors depend sensitively on the details of the compactification and we 
leave a careful determination to future studies. Indeed, the requirement T r < T max may 



significantly constraint the parameter space of the model [44]. 



5. Transfer of Energy into the SM Sector: D7 Wrapping a Non-Inflationary 
4-Cycle 

As discussed previously, the wrapped D7 scenario of the previous section may lead to 
dangerous ^-corrections to the inflaton potential. Thus, it may be desirable to exclude such 
a wrapping. Doing so, of course, also forbids a direct coupling between the inflaton and the 
SM fields. In this case the inflaton may still decay to SM particles via some intermediate 
bulk states. Such decays which involve bulk modes must necessarily be suppressed by 
the compactification volume, which is exponentially large. Here we discuss one possible 
scenario of this type. 

We could imagine the SM living on a D7 which wraps the big cycle T\, however, 



this leads to unnaturally small gauge couplings [37| and hence is phenomenologically dis- 
favoured. Instead, let us suppose the SM lives on a D7 wrapping some stabilized non- 
inflationary cycle, T3, which is stabilized to a value (T3) = 0(1). This wrapping is illus- 
trated schematically in Fig. [l] as "scenario 2" . 

At this point we must extend the discussion from section [4.3| to the case with more 
than two moduli fields, namely, instead of 2 fields T& (= t\) and r s (= T2) (corresponding 
to canonical moduli \ and f ) we consider 3 fields ri, T2 and T3 (corresponding to canonical 
moduli X; ^2 and $3). 

As we have seen, the endpoint of inflation is marked by the violent production of 
fluctuations 5t2 of the inflationary 4-cycle, T2. Near the minimum of the potential these 
T2 fluctuations mix with the total volume t\ and the non-inflationary hole T3 through an 
off-diagonal mass matrix similar to (|Oj) , except now there are three relevant moduli rather 
than two. Diagonalizing the mass matrix one finds three canonical fields x an d &i {i = 2, 3) 
which are schematically related to the 4-cycle volumes as 

&n ~ 0(V 2 / 3 ) x + 0{ y l/& ) ®i 

i 

6n ~ 0(V 1/2 ) + 0(1) X + Yl 0(V~ 1/a )*j where i = 2, 3 (5.2) 
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so that x is mostly the total volume 5t\ while &2, $3 are mostly the blow-up modes 8T2 
and 5t3 respectively. The modulus x is light, having mass 

m * ~ V3/2(in V )i/2 Mp (5 ' 3) 
while the moduli &2, $3 are heavier 

InV 



V 



M p (5.4) 



These moduli couple to gauge fields living on the D7 wrapping T3 via the interaction 
r^F^F^ , which gives us 

Ant = ^f 1 X F^F, V + Yl *< F ^ (5-5) 

i=2,3 

using equations (|4.5| ) and (|5.2|). Physically, these couplings originate from the mixing 
between T2 and t%, T3. From equation (5.2) we can estimate the magnitude of the moduli 
couplings. The largest coupling is, obviously, the one involving $3 since the SM D7 wraps 
the cycle T3. This coupling is set by the string scale 

1 V 1 ! 2 . . 

A ^ ~ ~rn s ~ ~Mp~ (5 - 6) 



A ^ ~ mHuV (5 - 7) 



On the other hand, the coupling to the large cycle is 

1 

x p 

where the factor M~ l comes from that fact that x is a bulk mode and the factor of (In V) -1 
is nontrivial [37|. The coupling between the inflationary 4-cycle T2 and brane-bound SM 
gauge fields is 

11 

A * 277 ~i^vv^ (5 ' 8) 

This is even more than Planck suppressed. 

Thus, the picture of reheating in this scenario is the following. Preheating produces 
copious amounts of 5t2 fluctuations which lead to a large number of canonical $2 moduli 
near the minimum, plus a small admixture of light x an d heavy $3 moduli. The $3 particles 
are the first to decay, they produce brane-bound SM photons via the process <3?3 — * 77 with 
a rate identical to (4.13): 

(InV) 3 



r$ 3 ^ 77 -s— M p (5.9) 



This decay, however, does not correspond to true reheating since at this point the energy 
density of the universe is dominated by the almost nonrelativistic 9 $2 particles with matter 
dominated equation of state, rather than SM radiation. Thus, the SM radiation which is 



See subsection 3.2 for a discussion of the fact that the bulk of the inflation fluctuations produced by 
preheating are nonrelativistic near the minimum. 
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produced at by the decay of $3 is rapidly diluted, p va dn ~ o -4 ; relative to the nonrelativistic 
$2 particles whose energy density dilutes as: p$ 2 ~ a -3 . 

True reheating occurs when the $2 particles subsequently decay into SM states via 
the suppressed coupling (|5.8[). This decay proceeds with rate: 



(lnV) ; 



It is the decay rate ( |5.10| ) which determines the reheat temperature of the universe: 

(InV) 3 / 2 



T r ~ O.l^TMp ~ 0.1-^ — ^ — M p (5.11) 

This is smaller than the result of the previous section, equation ( f4.20| ), by a factor of V -1 
corresponding to the fact that the $2 is a bulk mode in this set-up with suppressed coupling 
to brane-bound states (|5.8|) . This suppression is actually favorable since it keeps the reheat 
temperature ( |5.11| ) well below the decompactification scale T max , given by Q4.22| ). Taking 
V ~ 10 6 we have T r ~ 10 7 GeV and taking V ~ 10 9 we have T r ~ 10 GeV. Hence, for 
extremely large values of the compactification volume, V ~ 10 9 , the reheating temperature 
is so low that it may be difficult to realize baryogenesis at the electroweak phase transition. 
We leave this issue to future investigations. 

In this scenario the gravitino decays too early in the history of the universe to interfere 
with BBN: for V ~ 10 6 — 10 9 we have m 3 / 2 ~ 10 9 — 10 12 GeV. As previously, the residual 
volume modulus x particles decay after reheating but before BBN. 

The results of this section may also have some relevance for the scenario where the SM 
cycle is much smaller than the string scale, see [43]. We leave a detailed investigation to 
future studies. 



6. Stringy Reheating via Kahler Moduli Shrinking 

In this section, we propose yet another reheating mechanism in Kahler moduli inflation 
models. This mechanism can operate even in the absence of a D7 wrapping the inflationary 
4-cycle. Hence, in the ensuing text we suppose (as in section |5|) that the SM lives on a D7 
wrapping T3. 

The inflationary dynamics of r{t) corresponds to the shrinking of the 4-cycle associated 
with Ti . If this inflationary cycle of the internal CY shrinks to a minimal size r m comparable 
with the string scale, supergravity description breaks down and new, stringy degrees of 
freedom must kick in. A natural candidate for such degrees of freedom is winding modes: 
closed strings with a nonzero winding number with respect to the inflationary 4-cycle. 
Because of winding number conservation winding string are created in pairs with winding 
numbers of equal absolute value but having opposite signs. Two such closed strings can 
merge into a one with zero winding number, i.e. a string that can move away from the 
inflationary 4-cycle and into the bulk region of the CY. In other words, the energy of 
the shrinking inflationary cycle {i.e. inflaton kinetic energy) is transferred into closed 
string excitations living in the bulk. This is somewhat similar to the energy transfer from 
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annihilating branes into closed string excitations that occurs at the end of stringy warped 
throat brane inflation j|, |6| . The excited closed strings decay via many cascades into closed 
strings with the lowest level of excitation, i.e., KK gravitons. The produced KK gravitons 
interact with the SM brane scalar moduli Yi (corresponding to the transverse fluctuations 
of the brane) and also SM particles. Consider, for example, the interaction between KK 
gravitons and brane moduli Yi. This interaction is gravitational and can be described by 
a vertex jj T h*fJ < T fJ ' l/ (Y), the corresponding decay width is 



'p 



r^-^jf-. (6.D 

As discussed in || , this coupling is sufficiently generic to work for any low-lying degrees of 
freedom on the SM brane world. (In particular, world-volume fermions.) For the purposes 
of making a rough estimate of the reheat temperature we can safely assume that Tkk-^y ~ 
Fkk^sm- 

Let us now estimate (|6.1|). The KK modes have typical mass ttikk ~ 1/Rcy, where 
the overall size of the internal CY is Rcy ~ V l / & \[a' \ In the models under consideration 
we have 

i—. i Vv 



(6.2) 



m s Mp ' 

so that Rcy ~ 1F~> m KK ~ tH^j anci from ( |6.1| ) we have 

FkK->SM ~ -yf- • (6.3) 

The decay rate (|6.3j ) is suppressed as compared to the case with the wrapped D7, equation 
(4.13) by a logarithm which gives a numerical factor of about 1CP 3 . 



The decay rate (|6.3| ) allows us to estimate the reheat temperature as 

T r ~ 0.1^M p (6.4) 

Taking, for illustration, V ~ 10 6 we have T r ~ 10 11 GeV and taking V ~ 10 9 we have 
T r ~ 10 s GeV. As in the scenarios discussed in sections |] and || such a high reheat 
temperature is phenomenologically sensible because the gravitino is extremely massive. As 
in the scenario discussed in section || (but unlike the scenario is section |j) this reheat 
temperature is well below the decompactification scale, T r /T max < V _1//4 -C 1. 



7. Summary and Discussion 

Modular inflation, and other string theory inflation models, provide a natural playground 
for studying reheating in a context where (at least in theory) one can actually determine 
all decay channels of the inflaton to the visible (SM) sector from first principles. Here we 
investigated in detail reheating after modular/Roulette string theory inflation models in the 
context of large volume compactifications. Our results show that in realistic microscopic 
models the details of reheating can be rather complicated and may proceed through a 
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variety of channels, including perturbative decays, nonperturbative preheating and also 
inherently stringy processes. 

We found that modular inflation models (such as Kahler modulus or Roulette inflation) 
in particular display very rich post-inflationary dynamics. 

The initial stages of the decay of the inflaton proceed via extremely efficient non- 
perturbative particle production involving a combination of tachyonic instability and 
parametric resonance. This combination may lead to the most violent example of 
preheating known in the literature. 

The subsequent stages of reheating involve the transfer of energy from the inflaton into 
excitations of the visible sector. This phase is more model dependent since we must identify 
the location of the standard model of particle physics in the compactification volume. 
We have considered three separate scenarios for transferring the energy from the inflaton 
fluctuations St into SM degrees of freedom. These are as follows: 

D7 wrapping T2 (assuming a SUGRA description of reheating is valid, i.e., the min- 
imum hole size is larger than the string length): In this case the SM fields live on a 
D7 wrapping the inflationary 4-cycle, T%, so that there is a direct coupling between 
the SM fields and the inflaton r. Photon production (both perturbative and non- 
perturbative) dominates the decay of the inflaton. The resulting reheat temperature 
is given by fl4.2CP . Using V ~ 10 6 we obtain T r ~ 10 13 GeV and using V ~ 10 9 we 
have T r ~ 10 10 GeV. Although such reheat temperatures are extremely high, there is 
no problem with BBN because m 3 / 2 ~ 10 9 — 10 12 GeV is so large that the gravitinos 
decay before the onset of nucleosynthesis. However, such a high reheat tempera- 
ture may lead to decompactification. This scenario also has the disadvantage that 
^-corrections may spoil the flatness of the inflaton potential. 

D7 wrapping T3 (assuming a SUGRA description of reheating is valid, i.e., the mini- 
mum hole size is larger than the string length): In this case the SM fields live on at 
D7 wrapping some non-inflationary 4-cycle, T3 (we exclude the big cycle T\ since this 
would give unacceptably small gauge couplings). Near the minimum of the potential 
the produced inflaton fluctuations 5t2 mix with the fluctuations of the big 4-cycle 5t\ 
and also of the small non-inflationary 4-cycle 5t% through a non-diagonal mass matrix. 
The decay of the 5t2 fluctuations into brane-bound SM states via Planck suppressed 
operators leads to reheating. This scenario evades any dangerous ^-corrections to the 
inflaton potential and the fact that reheating proceeds through Planck suppressed op- 
erators suppresses the reheat temperature to well below the decompactification scale. 
Taking V ~ I0 6 we have T r ~ 10 7 GeV and taking V ~ 10 9 we have T r ~ 10 GeV. As 
above, there is no gravitino problem because m 3 / 2 ~ 10 9 — 10 12 GeV. 

Stringy reheating: In this case closed strings are produced when the inflationary 4- 
cycle size becomes of order the string length. These closed strings cascade into KK 
gravitons which can subsequently interact with brane-bound SM fields (wrapping 
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some non-inflationary 4-cycle) via Planck-suppressed operators. Taking, for illustra- 
tion, V ~ 10 6 we have T r ~ 10 11 GeV and taking V ~ 10 9 we have T r ~ 10 8 GeV. 
There is, again, no gravitino problem because is so large. Decompactification 
is evaded due to the Planck suppression of the interactions between bulk and brane 
modes. This scenario also evades potentially dangerous ^-corrections to the inflaton 
potential. 

For convenience we summarize the rates for the dominant decay channel of the inflaton in 
the various scenarios in the following table: 



Summary of Dominant Decay Rates 


t/m p 


Scenario 


(InV) 3 
~ V*~ 


D7 wrapping T2 


(InV) 3 


D7 wrapping T3 


1 


stringy reheating from shrinking hole 



In all cases considered there is the disadvantage that SUSY is broken at too high a 
scale to explain the small value of the Higgs mass. This happens because the typical 
compactification volumes which are favorable for inflation V ~ 10 6 — 10 9 are much smaller 
than the value V ~ 10 15 that is favored for particle phenomenology [37|. See, however, 
1 43] for a scenario where V ~ 10 6 — 10 7 may be compatible with TeV soft terms. 

Our first scenario, the SM D7 wrapping T2, is afflicted by two potential complications: 
^-corrections may destroy the flatness of the inflaton potential and the reheat tempera- 
ture is so high that it may lead to decompactification. Both of these issues will require 
further study, however, at first glance this scenario seems disfavoured. On the other hand, 
the scenario with the SM D7 wrapping T3 evades both of these potential difficulties and 
therefore seem most promising. For such models there are again a number of directions 
open for further study. Also, it would be interesting to perform a more detailed study 
the topology-changing transition in the case where the inflationary hole size shrinks to the 
string length at the end of inflation. We leave these, and other, interesting issues to future 
investigations. 
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